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ABSTRACT 


The production of electron-positron pairs by single photons In magnetic 
12 

fields >10 G has been Investigated In detail for photon energies near 

2 

threshold ( -Wu > 2 me ) as well as for the asymptotic limit of high photon 
energy. The exact attenuation coefficient, which Is derived and then 
evaluated numerically, Is strongly Influenced by the discrete energy states of 
the electron and positron. Near threshold. It exhibits a "sawtooth" pattern 
83 a function of photon energy, and Its value Is significantly below that 
predicted by the asymptotic expression for the attenuation coefficient. The 
energy distributions of the created pair are computed numerically near 
threshold and analytic expressions are derived in the asymptotic limit. These 
results indicate that as field strength and photon energy Increase, it becomes 
Increasingly probable for the pair to divide the photon energy unequally. 

This effect, as well as the threshold behavior of the attenuation coefficient, 
could have Important consequences for pulsar models . 
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I . INTRODUCTION 


12 

The possible presence of superstrong nagnetic fields (> 10 G) around 
astrophyslcal objects such as neutron stars has TOtlvated extensive study of 
the physical processes which can occur In these environments. In particular, 
quantum-electrodynamical processes such as magnetic pair production and 
synchrotron radiation In the high field limit have come to play such a central 
role In pulsar models that a complete understanding of these processes Is 
vital to further theoretical development. In addition, unexplored 
applications of high magnetic field physics to other objects, e.g. gamma ray 
burst sources, require further Investigation of single photon annihilation and 
pair production. 

The creation of electron-positron pairs by single photons in strong 

magnetic fields has been studied In a number of Instances. The first detailed 

treatments were carried out independently by Toll (1952) and Klepikov (1954) 

who calculated the attenuation coefficient for photon pair conversion using 

the exact Dirac equation for electrons and positrons In a constant, uniform 

magnetic field, treating the radiation field of the photon as a first order 

perturbation. Their results, which were In mutual agreement. Indicated that 

even for photon energies well in excess of threshold (hio ■ 2mc ) extremely 

high magnetic fields (B > 10^ Gauss) were required for a non-negllglble rate 

of pair production. In the limit of high photon energies the rate was found 

2 

to depend on the parameter x = ('ttu/2mc ) (B sin where 0 is the angle 

2 3 

between the field and the photon wave vector, and B^^ -me /e ^ - 
13 

4.414 X 10 Gauss. The constant B^^ Is the critical field strength. In which 
the gyroenergy ^ij)^ of an electron (or positron) Is equal to Its rest mass. 
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Since changes In the electron's energy perpendicular to the field occur In 
quantised Increments of AE ~ sets a scale on which to ^uge tht 

Importance of quantum effects In magnetic fields* The calculated conversion 
rates Indicate that magnetic pair production Is not expected to occur with 

significant probability unless x ^ even though it is kinematically 

2 

possible for any photon with hu > 2mc to produce a pair, no iKttter how 
low X • 

These results were later reconfirmed and extended by Erber (1966), Baler 

and Katkov (1968), Rassbach (1971) and Tsai and Erber (1974). These authors 

were primarily concerned with the derivation of asymptotic expressions for the 

total pair production rates in the limit of low field strengths (B « B^.^) and 

photon energies much higher than threshold. However, at energies near 

threshold, quantum eifects due to the discreteness of the electron and 

positron energies in the magnetic field stvongly affect the properties of pair 

production. When the photon energy and propagation direction are such 

2 

that -hti) sin 0 > 2mc (a situation which can occur in pulsar cascades) the 
asymptotic limits become unreliable. Unfortunately, the properties of 
threshold magnetic pair production have not been sufficiently explored in the 
published literature and are discussed only qualitatively if at all (Toll 
1952). Furthermore, no quantitative treatment of the energy distribution cf 
the created pairs exists, even in the asymptotic limit. Toll (1952) and Arons 
and Scharlemann (1979) do give qualitative discussions of the differential 
conversion rates, concluding that the pairs should most probably share the 
parent photon energy equally, but do not derive expressions for the energy 
distribution . 

This paper is an Investigation of these several important areas where 
gaps exist in the literature on magnetic pair production and which are 
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relevant to aatrophyalcal applications: properties of pair production near 

2 

threshold, > 2qc , and the energy distribution of the pairs, both near 
threshold and In the asymptotic limit* We also review the previous results 
with an emphasis on their physical significance and attempt to define more 
clearly the approximations under which the asymptotic results have been 
derived. In Section II, we begin with a discussion of the kinematics of pair 
production idiich, aside from being of Interest In lt:.elf, greatly helps to 
clarify the meaning of the results which follow. The derivation of the exact 
pair production rate for both photon polarizations Is outlined In Section III, 
In which we also present numerical calculations of the rates and pair energy 
distributions near threshold. Section IV reviews the derivation of the rate 
in the asymptotic limit of large photon energy. We also give a comparison of 
the exact and asymptotic rates and derive the asymptotic pair energy 
distributions. In Section V, ve show how the results can be generalized to 
arbitrary photon directions and non-vanlshlng electric fields through Lorentz 
transformations and discuss in Section VI the importance of some of the new 
features for models of astrophyslcal sources. 

II. KINEMATICS OF MAGNETIC PAIR PRODUCTION 


Consider a photon with energy u> propagating at an angle 0 to a uniform 
magnetic field (Natural units, with +I ■ c • 1, will be used throughout). 
Without loss of generality. It Is possible to choose a frame In which Jl ■ 

A 

and k ■ u) (0, sin 0, cos 0). Then the energy-momentum conservation 
laws governing magnetic pair production may be written as 


(i) *■ E , + E, 
J k 


(la) 


u) cos 0 ■ p + q. 


(lb) 
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where the energy quantum numbers (j, k 0, 1, 2, •..) and the longitudinal 
laomentum components (p, q) refer to the positron and electron respectively* 

The Landau energy levels frr the two particles have the forms (see for example 
Johnson & Llppman 1949) 


r 2 2 - 1/2 

E - fp + m (1 + 2J -2—)] ' 

J B 

cr 


r 2 2 B , 1/2 

E. - fq + m (1 + 2k -2-)] ' 

B 

cr 


<2a) 

(2b) 


Note that the conservation laws do not require conservation of momentum 
perpendicular to the field direction. Indeed, the fact that the external 
field participates In the particle momentum transfer Is what makes this 
first-order transition kinematically possible. However, it Is the "rigidity” 
of the Intense field lines which preserves the energy conservation law [Eqn 
(la)]. 

In the derivation of the photon attenuation coefficients. It Is 
convenient to specialize further to the frame In which ■ 0, so that 
p ■ -q . In this frame, equation (la) may be solved for p to yield 


r 2 2 at: 1/2 

p - p(j,k) - ± m - 1 - (J + k) B» + (j-k) -2-] (3) 

2 

4o)’ 

9 

where u)' - w/2m, B' ■ • Since p must be nonnegative, the allowed 

transitions to e'*'e' pair states represented by Integer pairs (j,k) are just 
those points in the J-k plane enclosed by the curve p^ ■ 0, as shown In 
Figure 1. 
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It turnd out that tranaltlons to states close to this boundary are the 
most probable, as might be expected from the fact that classically the pair 
should most probably appear with high "forward” momenta (i.e., large 
transverse energy levels and relatively small momentum components along the 
field direction)* In any case, the maximum Landau level obtainable by either 

member of the pair is found by assuming that either J or k - 0 on the boundary 

2 

curve p (J,k) > 0, from which it is found that 


J 


max 


k 

max 


B' 


( 4 ) 


2 

For u)' » 1, this reduces to ~ C = 2u' /B’ . The parameter 5 appears 

in the derivation of the attenuation coefficient in Section III below. Similarly, 

2 

along the symmetry axis J - k, it is found that < (o)’ - l)/2n* . 

The number of Integer (j k) pairs enclosed within the boundary p^ ■ 0 Is 

2 2 

given approximately by (w*, B’ ) ~ 2u)* (oj’ + 2)(w' - 1) /3B' , which 

4 2 

for the limit o)'»l, reduces to (w',B') ~ 2u)' /3B' . Hence for 

fixed B*, an Increase in the photon energy dw' corresponds to an increment 
3 2 

dNgtates ~ duj'/3B' . This number Is just twice the number of 

singularities which appear in the photon attenuation coefficient between 


energies u' and u)' -i du', as will be discussed in Section III below. 


III. DERIVATION OF THE PHOTON ATTENUATION COEFFICIENT 


The S-matrlx element for magnetic pair production is simply the first 
order transition 


fd 


X 
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where the electron wavef unctions are solutions of the Dirac Haailconian for 
motion in a constant, uniform magnetic field* The precise form of these 
wavefunctions depends both on the choices of gauge and constants of the 
motion. The wavefuctions used in the present analysis are described 1 ' detail 
by Daugherty and Bussard (1980); for alternative solutions, see for example 
Johnson and Lippmann (1949) . 

The attenuation coefficients for photons whose electric vectors are 
polarized parallel or perpendicular to the magnetic field direction may be 
written as 


‘l.i 


- I I 

T j>o k>o 


!i iEft HjiM- h 

2it 2ir 2it X 


db 

~2 

2itX 


fi 


I. .1 


(5) 


where L T Is the space-time volume element and X ■* 1/eB. The summations over 
J and k range over all kinematically allowed final states (cf. Fig. 1). The 
s_ and 8^ summations range over the spin states of the electron and 
positron. The continuous variables (a,b) are eigenvalues of the x- coordinates 
of the orbit centers, which are constants of the motion for the trave functions 
used here. As indicated In the earlier discussion on kinematics, the sum over 

Landau levels j and k is restricted to those states for which 

2 

P (J.k) y 0 [see equation (3)]. 

The derivation of the attenuation coefficients using the wavefunctions 
specified above is straightforward and confirms the results urlglnslly 
obtained by Toll (1952) and Klepikov (1954) . For details on techniques for 
handling the matrix elements, see also Daugherty and Bussard (1980), Appendix I. 
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The results may be written In the foras 

{(E.E. -hi^SdMCJ.kif+lM (J-l.k-l)l^ 
2C J klpj^l ^ ^ 


2 

+ 2/W B’m [M+(j,k) M(J-l,k-l) + M+(J-l,k-l) M(j,k)]} 


(6s) 


R.(u’,B‘) - xf-T T|^|{(E E ■hIl%^(|M(J-l,k) I +|M(J,k-l)|^) 
’ j k'Pjk' j k 



(6b) 


where F, 


2 

2a)' 

B’ 


M(J,k) 


<-l) 


G-S 




(O, 


( 7 ) 


and G - max (j,k), S - min (j,k).. 

G”S 

L g (O l8 the generalized Laguerre polymomlal, defined according to the 
standard in Abramowltz & Stegun (1964) . These results may also be extended to 
the cases In which j or k (or both) Is zero, If one adopts the convention 
(O = 0. 

The attenuation coefficients are immediately seen to have singularities 
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or absorption edges whenever u* end B' have values such that the longitudinal 
mooentuB p(Jtk) vanishes for any of the allowed transitions (J,k) In the 
sunmatlons of Equation (6). If the attenuation coefficients are considered as 
functions of ci)' for fixed B', these singularities occur at the energies 

Jk "V2(*'1+2JB’ + /l+2kB') (8) 

Figure 2 shows this behavior In the parallel and perpendicular attenuation 
coefficients In the vicinity of the threshold energy for several field 
strengths. IWo oajor trends with decreasing field strength to be noted In 
this figure are the decrease of the typical spacing between peaks and the 
sharp Increase of the attenuation coefficient (averaged over saall energy 
Intervals) at energies near threshold. The spacing between peaks can be 
understood from the klnematical discussion In Section II, and In particular 
for u* » 1, the number of absorption peaks between a>' and u’ -f do)' Is 
approximately 


dN 


peaks 


- Va dN 

0 


to’ 


tates 


du' 

3 B’ 


The "sawtooth" behavior of the attenuatl'^n coefficients has been 
discussed both by Toll (1952) und Klepikov (1954), but has generally been 
ignored in astrophyslcal applications to date. One argument for neglecting 
such effects is that In pulsar magnetospheres one Is typically dealing with 
energetic photons (w' » 1) and field strengths somewhat less than B^,^. 

Hence one might expect sufficiently small spacing between peaks so that the 
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well known aiyiaptotlc, energy-everaged forw of the attenuation coefficlenta 
(sec Section IV below) are suitable approximations. However, it oust be noted 
that in these same situations, one typically finds that these photons 
propagate at small angles to the local field direction. Since the asymptotic 
forvs, like the discrete summations {cf. Eqn. (6)] above, apply only to 
ptopagation perpendicular to the field, it is in fact necessary to verify (by 
Lorentz transformations) whether it is still true that u* >> 1 in the frame 
for which This question will be pursued in Section V below, where it 

is shown that typically "energetic” gamma radiation may actually be subject to 
the near-threshold behavior of the attenuation coefficients. 

The relative individual contributions of the (j,k) states in the 
summations of Equation (6) are Illustrated in the plot of Figure 3. As noted 
earlier, the states nearest the boundary p^ - 0 are by far the most 
significant. However, Figure 3a suggests that there should be a broad energy 
distribution for each member of the created pair. This indication is 
confirmed by the energy distribution obtained by a sampling technique over the 
(J, k) plane, which is shown in Figure 4. For B - lO^^G, u - 3 >teV, there is 
a strong tendency for the pair to divide the phrcon energy equally. However, 
for the case B ■ a> ■ 3 MeV, there is actually a greater probability at 

each member of the pair will be created with an energy near either the minimum 
(E “ me ) or maximum (E ■ u) - me*) allowable. This behavior near threshold 
will be compared to energy distributions in two distinct asymptotic regions in 


Section IV. 
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IV . ASYMPTOTIC EXPRESSIONS 


a) Attenuation Coefflclenta 

For photon energlea an'l field strengtha auch that the parataeter 

2 

C “ 2b)' /B' (and hence the nuaber of allomd J-k atatea for the created pair) 
becodten large , it la poa^ible to reduce the attenuation coefficienta (Eqn (6)' 
to aaymptotic foroa which are effectively averaged over aiaall energy intervala 
to remove the sawtooth behavior. This la feuslMe because the number of peaks 
in such small intervals increases rapidly beyond threshold, while the peaks 
themselves become more narrow at an even taster rate (Toll 1952) . The 
averaging procesii is actually done by replacing the discrete sumations over 
j,k (in which individual singular contributions make the entire sum infinitti) 
with Integrations over continuous variables. Since the contributions from 
transitions on opposite -<dea of the symmetry axis J • k in Figures I and 3 
are identical, it is reasonable to start by defining such variables as 
follows: 

u - (J + k)/r, (9) 

V - (J - k)/r. 

In t:^rms of these quantities, 

2 2 , 2 2 , 

p ■ m f b)' (l-2u+v )-l 1 (lO) 

Since one Integration boundary Is determined by the curve p‘ - 0, It proves 
convenient to define a further variable (cf. Figure 1) 
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♦ 


2 1/2 

(1-2U+V ) ' 


( 11 ) 


2 2 2 2 

which Is then related to p by p ■ ■ [u* ♦ " 11 • The wost probable final 

states lie in the vicinity of ^ ■ 1/m' , hence it will be reasonable to 

t&in 

expand quantities in the Integrand In powers of 

The principal ingredient in obtaining these asymptotic results la a 
suitable approximation for the matrix element M(j,k) defined in Equation 
(7). A derivation of asymptotic forms for M(j, k) based on that of Klepikov 

(1954) la given in the Appendix. The reduction of the attenuation 

coefficients to their as)TBptotic forms than begins by rewriting the matrix 

elements M (J, k), M (j-l, k-l), M (J-1, k) and M (J, k-1) entirely in terms 

of ^ and V, using Equations (A3), (A6) and (A*^) . Also, with the use of 
Equations (9) and (10), the various energy terms in Equation (6) may be 
expanded to give the following: 

2 2 2 , 2 2 2 , 

EjEj^ m - p ~ m f2+u* (1-v - <t> ) I (12a) 

2 2 2 2 2 2 

EjEj^ + m + p ~ m ( 1 )' (1-v + <j> ) (12b) 

^ 2 2 2, 2 2 n + , 

2/jic B* a ~ m u' fl- V - ] (12c) 

(1 - vS 

The insertion of these expressions and those for the matrix elements into 
Equation (6), an. '.le replacement of the suras by Integrals, 

1 l-(|> 

2 

y y -*■ 2 f d<^ r av (zc ♦) 

J k ^ o 

M* 
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leads to expressions for the (sooothed-out) attenuation coefficients valid for 

Replacement of the variable ♦ by n - allow® the attenuation 

coefficients to be written In the forsa 


, , 1° _ « — r d„ f 


' * , 2 .2,, 1 /n'-l o J-v2 

3w (1) B 


2 4 


’{fn + n — 2 ^2/3(5)^ 

2 ' 1- V 

(1- V ) 


(13a) 


p o 8 


n r dv 


fdrv-^ f 


-«2 


* 3,2 ^ B- > ° 2-v 


** 2 4^2 

•| — ^ Ki ^3 (O + n — j K 2/3(C)l 


_n 

1-v 


(1-V ) 


(13b) 


where ^ 


3 



1-v 


Klepikov (1954) has written the attenuation coefficient for unpolarlzed 
photons In the form 


- (R, ^ Rp - 

2 ^ 


o B' --- 


2 2 
3it X 


f“dx f*dyf2 cosh y cosh^ x 


o 


o 
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2 2 2 3 2 322 2 3 

Kl i 3 i — co*h y co#h x) - slnh x co«h x Ki# 3 ( — coah y coah x) 

' 3x ' 3x 

2 5 2 2 2 3 , 

+(2 coah y- 1 ) cosh x K 2/3 ( — cosh y cosh x) | (lA) 

3x 

where x ~ . This expression Is found to be equivalent to the average of 

equations (13a) and (13b) through the replacements 

cosh X ■ n 

cosh y ■ — L— . ( 15 ) 

2 lu 

In these variables, sinh x"(n-l) ^ ■ p/m, so that the x integration in 
Equation (14) is essentially an Integration over the longitudinal momentum. 

The variables v in Equation (13) and y in Equation (14) are related to the 
energy of the electron or positron through Equations (2), (9), (10) and (11), 
which give: 

V “ 1 1 - 2 e I 
2 

1 - V - 4e(l - e) (16) 

where t - Ej/w Is the fractional energy received by one member of the 

pair. Therefore, the v and y Integrations in Equations (13) and (14) are 
Integrations over the pair energy. 


Equations (13) may be further reduced In the limits x « 1 »nd x 1* 
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For the case of soall x particular, both of the modified Bessel functions 
approach the asymptotic form 

*1/3 (C) ~ *2/3 (?) ~ «■' • 

Furthermore, the Integrand becomes sharply peaked around the symmetry axis 

2 2 

(v-0), so that It Is permissible to approximate 1/1-v as 1+v in the 
exponential above. Integration over v after this step, followed In turn by 
the approximation n ~ 1+p, P<<1, yields the well-kno»m result 

a A 

R ~ 0.23 B' exp (- — ) . (17) 

* 3x 

This is In fact the expression usually assumed in astrophysical calculations, 

although as shown in Section V below It Is sometimes Inappropriate. The 

convergence of the exact attenuation coefficient derived in Section III (cf. 

Eqn. (6)] to this asymptotic limit Is Illustrated In Figures 5 and 6. The 

exact attenuation coefficient is at least several orders of magnitude less 

^12 

than the asymptotic value near threshold for field strengths around 10 
Gauss. Below threshold, of course, the asymptotic expression becomes 
completely invalid. 

b) Pair Energy Distributions 

The energy distribution of each member of the created pair (which are 
identical because of the symmetry of the Integrand about v ■ 0) can then be 
obtained by Integrating Equation (13) first over n, and then substituting for 
V In terms of e in the resulting integrand. Tsai and Erber (1974) have shown 
that the polarization-averaged attenuation coefficient [Equation (14)] can be 
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expressed in a sinpler form, 

R , fH f dv K 2/3fV3x(l - vS 1 . (18) 

^ ISiTx o ^ ^ 

By the simple change of variables of Equation (16) we obtain the differential 
attenuation coefficient as a function of pair energy: 


dR(c,x) 

dc 


V /~T~ 12 + e(l-e)] ^ 

ir 9irx e(l - e) ^/3 


tl/3xe(l-c)l 


(19) 


which has the asymptotic forms 


"• ot B ^ 

^ 9. L_tL±^L\l£>l_ exp [-l/3xe(l-E)], x « 1 (20a) 

^ 3/~ Jirx *' eTl-e) 


“8' / ot 

o .426 


1/3 




[2 + c(l-E)] 
1/3 

[e(l-e)] ' 


X 


» 


1 

€(l-e) 


(20b) 


These energy distributions have a number of interesting properties. The 
first to note is that the parameter x alone determines the shape of the 
distribution, which has been plotted in Figure 7 for different values of x 
using the exact form in Equation (19). For small x . the distributions are 
centrally peaked around e • V 2 and have half-widths ~ .36 / x X 

For large x » the distributions completely change shape, actually becoming 
peaked at e - 0 and e “ 1. For x » 1« the shape of the distribution 
becomes independent of x altogether. Near e •• 0 and e ■ 1, the 
condition X 1 "°t Insure that the argument of K 2 / 3 (x) is small, and 

the exact expression of Equation (19) is needed to describe the drop to zero 
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at the edgea of the distribution. By comparing the curves In Figure 7 to 
those of Figure 4, it Is apparent that the energy distributions near threshold 
have a similar dependence on x • 

From the behavior of the pair energy distributions, «e can conclude that 
In the case of low fields and photon energies near threshOi.d the pair tend to 
share the photon energy equally. However, In the case of large fields or high 
photon energies, there is an Increasing tendency for one member of the pair to 
take all the energy of the photon, while the other receives little more than 
its rest mass. In the limit x this Is actually the most probable final 

state of pair production. 

V. LORENTZ TRANSFORMATIONS OF THE ATTENUATION COEFFICIENTS 


The results of the previous sections all apply to the special Lorentr 
frame in which Jt.JS ■ 0. To generalize the (unpolarlzed) photon attenuation 
coefficient to frames In which the photon wave vector makes an arbitrary angle 
9 with the magnetic field, tt Is sufficient to perform a Lorentz 
transformation with velocity v ■ cos 9 along the field direction. The 
transformation law for the attenuation coefficient is most easily understood 
In terms of its Inverse, the mean free path, which obeys the simple Invariance 
law 


2 2 2 2 

- t - X - t (21) 

o o 

From this law and the transformation properties of w and B, it is immediately 
possible to write 
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ii (b), B) - tin e R (u> sin 6, B) (22) 

- o “ 

where R Is now used to denote the functional fora of the attenuation 
o 

coefficient in the special frame for which ■ 0, as derived in previous 
sections. Note that this transformation law is valid for both the exact forms 
in Eqn (6) and the asymptotic limits of Section IV. 

In partic'ilar, for the asymptotic case x th® well-known 

approximation [Equation (17)] for the attenuation coefficient appears to be 
modified simply by the replacement of B with B sin 6. However, from Equation 
(22) it may be seen that it is really id which becomes id sin 6 in the 
product X “ (‘*)/2m)(B/B^^) . This distinction Is quite important, since in the 
non-asymptotlc case the attenuation coefficient does not depend only on x* 

The behavior of the transformed coefficient near threshold is illustrated 
in Figure 8. Note that the threshold energy itself is displaced upward by the 
factor 1 / sin 6, which also is the factor by which the energy scale is 
stretched in the new frame. The reduction of the attenuation coefficients by 
the factor sin 6 is also evident. 

The transformation property implies an additional point of Interest in 

2 

pulsar cascade theory. Here a primary electron with energy ymc , accelerated 
along a curved field line, may emit gamma radiation nearly along its direction 
of motion, within a forward cone of half-angle 0 ~ 1 /y. To follow ouch a 
photon along its trajectory outward through the magnetosphere, it is necessary 
to perform a sequence of Lorentx transformations at small intervals, where the 
field makes a slowly Increasing angle with the photon wave vector. As Figure 
8 illustrates, this sequence of transformations may initially put even high- 
energy photons below the threshold for pair production, and in such cases 
these photons will propagate at least until they exceed the local threshold 
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energy. Hence the near-threehold forme of the attenuation coefficients must 

be used over some portion of the photon trajectory, and the commonly used 

asymptotic form may not be appropriate over any portion of the actual photon 

trajectory before Its annihilation. 

It may also be noted that a Lorentz transformation In which ^ la 

perpendicular to ^ may also be used to generalize the attenuation coefficient 

to situations In which both electric and magnetic fields are present, at least 

2 2 

within the restrictions B-E>0, B»E"0 (Daugherty and Lerche 1975). 

Such cases are Important for tracing cascades In rapidly rotating pulsars (as 
for example the Crab pulsar) In which significant rotation-induced electric 
fields are expected. Here the transformation law for the attenuation 
coefficients may be written In the form 

5(o), B) - y( 1 - Ti^v) R^(yu(1-ti^v), ®) (23) 


in which (n^, Tiy, n^) are the photon direction cosines in a frame chosen so 

A A 2 

that B»Bz, E«Ey, the velocity of transformation y ■ bxB/B , and the 

2 2 _ 1 /, 

Lorentz factor ^“(l-E/B) 

The full generalization to situations in which electric fields both 
parallel and perpendicular to ^ are present requires explicit calculation of 
the attenuation coefficients using the Schwinger proper-time technique 
(Urrutla 1978). However, for known pulsars the ratio E^/B is sufficiently 
small that such cases are of only minor Interest in cascade theory (Ayasll 
1978, Harding et al. 1978). 
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VI. DISCUSSION 

We have obtained results on several aspects of magnetic pair production 

which should be of Interest In future astrophyslcal applications. The pair 

production attenuation coefficient near threshold Is dominated by magnetic 

field quantum effects, showing sharp resonance peaks at photon energies 

corresponding to discrete states of the electron and positron. As the field 

approaches the critical field strength, the spacing of these peaks at 

threshold Increases, approaching me /sln6. As the photon energy Increases 

above threshold at a given magnetic field strength, the peak spaclngs decrease 

rapidly and the mean value of the attenuation coefficient approaches, from 

below, the asymptotic limit derived for very large pair quantum numbers. The 
2 

parameter C = 2u’ /B* , which Is the number of the highest energetically 
allowed Landau state available to either the electron or positron. Is not a 
Lo rents Invariant quantity and must be evaluated In the frame In which the 
photon propagates perpendicular to the field. Thus, the asymptotic 
condition ^ »1 does not ensure that threshold (quantum) effects are 
unimportant In an arbitrary frame . Neither Is the threshold energy for pair 
production frame Invariant; In an arbitrary fratne in which the photon 
propagates at an angle 9 to the field, the threshold energy Is 2mc /sln6. 

This result has the somewhat surprising consequence that even extremely high 
energy photons. If they are propagating at very small angles to the field, may 
have perpendicular energies below threshold. 

Since the complicated expressions [cf . Eqn (6)] for the threshold 
attenuation coefficients are unwieldy for practical calculations Involving 
pair production, we give here an analytic expression which approximates their 
average behavior, and also describes their approach to the asymptotic limit. 
For x< 1, 


pAdC 

quality 
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R « 0.23 |f B’ exp f- - V . > ] m*> i, 

0 u»'< 1, 

-2 .7 - .0038 

where f(w'.B’) - I + 0.A2 id’ B* (24) 

The function f(b>',B*) has been obtained by a three paraaeter fit to the 

attenuation coefficient R((i>',B') averaged over fixed Intervals of o)' large 

enough to contain at least several resonance peaks. The fit Is not terribly 

sensitive to the B* dependence, owing to uncertainties from averaging over the 

the peaks. We point out that the threshold condition, u’H, must be put In 

artlflcally, and that the expression for arbitrary photon angles may be 

obtained using eqn (22) . It Is Interesting tc note that while the asymptotic 

expression (Eqn (17)] Is a function only of x» exact attenuation 

coefficient, as approximated by Eqn (24), depends on u* and B' Independently. 

We have also found that the probability distribution of pair energies Is 

determined by the parameter x> which ^ Lorentz Invariant (for the case of a 

pure B field and transformations along B^) . For small Xt probable 

pair energy Is half of the parent photon energy. In the limit of large x> 

most probable final state Is one where the photon energy Is given almost 

entirely to one member of the pair. This unequal division of the photon 

energy was shown by Sokolov et. al. (1974) to be the case for B > B^^ where, 

of course, x 1 b always greater than 1. From the results derived in Section 

IV, however, we have shown that unequal dlvle .on of the photon energy can 

occur also for B < B , as long as x ^ !• Probability distributions of pair 

cr 

energy very similar to those of Figure 7 have been computed for both y-y pair 
production (Burns and Lovelace 1982; Burns, private communication) and y- 
nucleus pair production (Rossi and Greisen 1941) which also are seen to change 
their shape as a function of a single parameter characterizing each process. 
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The production of pairs with one member of the pair receiving all the 
photon energy, and therefore occupying a high Landau state, i^ile the other 
member receives only its rest mass, thus occupying the ground state, is 
possibly related to the dominance of ground state transitions in synchrotron 
radiation* It has been found (Sokolov et al» 1973; White 1974, 1976) that 
transitions of an electron from an excited state to the ground state in a 
magnetic field become more probable as the field approaches and exceeds the 
critical field strength. This tendency of the electron to radiate away all 
its energy via a single photon causes the tall of the synchrotron spectrum to 
rise near the "tip", where flui equals the electron energy. Although the 
kinematics of this process are different from those of pair production, the 
probability distribution of initial and final states of synchrotron radiating 
electrons in a near critical field would look very similar to the pair 
distributions we have found for \ » 1 . 

The properties of magnetic pair production explored in this paper have 
immediate significance for pulsar models. In current polar cap models (eg. 
Ruderraan and Sutherland 1975, Arons and Scharleraann 1979), particles are 
accelerated to high energies (lO^^-lO^^ «V) along magnetic field lines and 
emit Y“ray photons through curvature radiation. These photons are initially 
propagating parallel to the field and thus have w sln0"O. In order to produce 
pairs, they must acquire a sln0 large enough to (at least) reach threshold by 
traveling a straight trajectory in the curved and rotating dipole field. If B 
^ 10^^ G, these photons will have o) sin0 ~ 2mc^ when they pair produce and the 
asymptotic attenuation coefficient will under-estimate their mean free paths 
because (1) it overestimates the exact attenuation coefficient near threshold 
and (2) it has a non-zero (and possibly large) value below threshold. The 
absorption of photons In strong fields may be further complicated by threshold 
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behavior of the vacuum index of refraction. Shabad and Uaov (1982) have 
argued that curvature photona in fielda > 5 x 10^^ G will be refracted along 
the direction of B. If thia ia the caae, then the photon mean free patha will 
be even larger. In any event, the increaae in mean free patha from the pair 
production threshold behavior may be regarded aa a lower limit. 

Since the e'*‘e~ pairs produced within the accelerating region tend to short 
out the electric field, the spatial extent of this region is proportional to 
the photon mean free path. Larger mean free patha will increase the distance 
over which the voltage drop occurs. However, due to the weak dependence of 
the size of the voltage drop on the altitude at which the shorting-out occurs 
in these models, the acceleration energies may not be much larger. 

The increased mean free paths are expected to have a larger effect on the 
cascades which occur above the accelerating region. Threshold effects on 
individual mean free paths will be magnified by the numbers of photons 
involved and could produce significant changes in the predicted y-ray and pair 
spectra for fields in excess of ~ 4 x 10^^ G. Use of the computed pair energy 
distributions In simulations of pulsar cascades may also have interesting 
effects. Most previous calculations of this sort (Daugherty and Harding 1982, 
Ogelman et al. 1976) have assumed that the pairs are always produced with one- 
half the parent photon energy. As we have shown here, the energy 
distributions are in fact quite broad for x ^ 0.1 when pair production ia 
most likely to occur. Both of these topics will be the subject of future 
investigation. 
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APPEHDIX 


Derivation of Asymptotic Forms for M(j,k) 


Following Klepikov (1954), we start froa the contour Integral 
representation for the generalized Laguerre polynomial 



1 

■STl 



C* (1-2) 

(-Z) 


G 

S ’ 


(Al) 


where the contour encloses the origin. Replacement of G and S by the 
paraaetlzed variables defined In Equations (9) and (11) above allows one to 
write the saddle points of the Integrand as 


2^ • V2 (1 - ± ♦) (A2) 

Both of these points lie on the positive real axis, and the Inner point z. may 
be crossed by Initially deforming the original contour about the origin. 

However, such a straightforward application of the method of steepest 
descents proves Insufficient when ^ ” 1/u', for then the two saddle 

points (which are oriented In mutually perpendicular fashion) im>ve quite close 
together . Hence the precise choice of contour through the point z_ Is 
critical. Klepikov's approach Is to make the transformation z “ 2 _e^ and 
then to expand the resulting exponential function of t about the origin (which 
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% 

is now also tho saddls point) • ratalning powsrs up to t • This asthod ylalds 
an approxioation for tha Laguarra polynoaials which, whan coablnad with 
Stirling *s approximation for tha factorial functions in Equation (7), yialds 
finally 

|M(j,k)| > !m(^,v)| l-ri.3(;) (A3) 

/3T /T^2 ' 


where 


C 



(A4) 


An unfortunate complication which arises in the expressions In Equation 
(6) for the attenuation coefficients is that the quantities M(J>1, k-1), M(J- 
l,k) and M(J,k-l) all appear in addition to M(J,k). While in the limits of 
large (j, k) these adjacent matrix elsMnts approach M(j,k) itself, the small 
remaining differences must be taken into account in Equation (6)« The 
principal reason is that each coefficient is found to be the difference of two 
terms which are nearly Identical, so that the small differences in the matrix 
elements in distinct terms are significant in the final result. The technique 
for binding the adjacent matrix elements involves the use of the relations 


dL® (x) 

X 2 nL®(x) - (n+a)L® .(x) 

dx ® ^ 


<A5a) 


V— - - L , (x) 

dx n-1 


(A3b) 
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l“"^(x) - L® (x) - L® , (x) 
n n n-1 


(A5c) 


for the generalized Laguerre polynoalala to prove the following relations 
bet%feen the matrix elements: 


/jk M(J-l.k-l) - 

- (J+k-^) M(J,k) -C 

(A6a) 


2 

dK 


i/jC M(j-l.k) - ■ 

- f^“^'*'^lH(J,k) - 1 

(A6b) 


2 

dK 


M(j,k-1) - 

(^■1|lL.)M(J.k) + C 

(A6c) 


From Equation (A3) It Is found that 


dM(J,k) 
dC “ 




2 

-L » 


K2/3(0 


1A7) 


Substitution of Equations (A7) and (A3) Into (A6) then gives expressions for 
the adjacent matrix elements in teriM of i and v. 
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Figure 1: 


Figure 2: 


Figure 3: 


FIGURE CAPTIONS 

Kinenatlcally allowed Landau states for the created electron-* 
positron pair with principal quantua numbers (j, k), for the 
case B “ u - 3 MeV. Ihese states are deterwlned by the 
condition p^ > 0 [cf. Equation (3)] which, as Illustrated, Is 
the area Inside the line p ■ 0. The variables v and A ere used 
In the derivation of the asymptotic limit In Section IV [cf. 
Equation (9)]. 

Exact attenuation coefficient for pair production [cf . Equation 
(6)] In the case where the photon propagates perpendicular to a 
pure magnetic field, plotted against photon energy. Threshold 
energy is 1.022 MeV, below which R(u) ■ 0. 

a) Attenuation coefficients for polarization (of the electric 
vector) parallel and perpendicular to a magnetic field at the 
critical field strength. 

b) Attenuation coefficient averaged over both polarizations 
for B - 5 X lO^^g, 

Relative contributions of Individual (j, k) states of the 
created electron-positron pair to the total pair production 
attenuation coefficient for the c.ises 

a) B » B^,j, ■ 4.414 x lO^^G, w - 3 MeV, having - 11 and 

b) B - lO^^G, 0 ) - 3 MeV, having » 502 

The vertical scales are log scales with an artificial lower 


cutoff . 
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Figure 4: 


Figure 5: 


Figure 6: 


Figure 7 : 


Energy distribution of one waber of the created pair, conputed 
fr<»B an Integration over the (j,k) probability distributions of 
Figure 3* Plotted here Is the normalized probability density 
as a function of pair energy E divided photon energy (d. Any 
asymaetrles In the distributions about E/u) >0.5 are due to 
statlsltcal errors In the Hante Carlo sampling used. 

Comparison of the exact and asymptotic pair production 
attenuation coefficients, R and R^, averaged over both photon 
polarizations, for two field strengths. Each "X" point for B ■ 
lO^^G Is an average of the exact attenuation coefficient over 
the many absorption edges (which are much more closely spaced 
then the absorption edges for B ■ 5 x 10 G) In the surrounding 
energy Interval. 

Convergence of the exact attenuation coefficient to Its 
asymptotic limit for a field of lO^^G. The ratio of the 
asymptotic to the exact attenuation coefficient, as shown In 
Figure 5, are plotted against photon energy. Dependence on the 
parameters x ^ (see text) Is also shown. 

Energy distribution of one member of the electron-positron 
pair, derived In the asymptotic limit [cf . Equation (19)] for 
different values of x “ ('ttu)/2mc ) (B/B^j.). The quantities 
plotted here are the same as those of FIG. 4. 


Figure 8: 


Lorentz-transformed pair production attenuation coefficient for 
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the case sln6 ■ 0.2, where 6 Is the angle between the magnetic 
field and the photon propagation vector. 
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